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Fluctuations of electric and magnetic fields in the coUisionless plasma are found as a solution of the 
initial value linearized problem. The plasma initial state is on average stationary and homogeneous. 
When the state is stable, the initial fluctuations decay exponentially and in the long time limit 
a stationary spectrum of fluctuations is established. For the equilibrium plasma it reproduces the 
spectrum obtained from the fluctuation-dissipation relation. Fluctuations in the unstable two-stream 
system are also discussed. 

I. INTRODUCTION 

Spectrum of electromagnetic fluctuations is an important plasma characteristics studied in various contexts. In 
terrestrial experiments the spectrum, which is observable through scattering measurements, signals, for example, an 
onset of plasma instability or turbulence. Electromagnetic fluctuations in primordial cosmological plasma are analyzed 
to explain an origin of magnetic fields in the Universe. 

The fluctuations can be theoretically described using several methods reviewed in the classical monographs P, 0] ■ 
Modern field-theory techniques developed for relativistic plasmas are worked out in [3(, ^ . Physically most appealing 
seems to us the method proposed by Rostoker [s] and Klimontovich and Silin which is clearly exposed in the 
handbook Q . The method, which is applicable to both equilibrium and nonequilibrium plasmas, provides the spectrum 
of fluctuations as a solution of the initial value (linearized) problem. The initial plasma state is assumed to be 
on average stationary and homogeneous. When the state is stable, the initial fluctuations are explicitly shown to 
exponentially decay and in the long time limit one finds a stationary spectrum of fiuctuations. In this way one obtains 
for the equilibrium plasma the spectrum which is alternatively provided by the fiuctuation-dissipation relation. When 
the initial state is unstable, the memory of initial fluctuations is not lost, as the unstable modes, which are present in 
the initial fluctuation spectrum, exponentially grow. 

The fluctuations of the distribution function, electric charge or longitudinal electric field can be found rather easily, 
see 0, 0| • Analytic computation of the magnetic field fluctuations appear to be quite lengthy and tedious while 
the computation of fluctuation spectrum of the electric field, which is not constrained to be purely longitudinal, is 
a real challenge even in the coUionless plasma, as one has to take into account and sort out numerous terms. Up 
to our knowledge such calculations have not been published. In this article we study the fiuctuation spectrum of 
electric and magnetic fields in detail. In the case of equilibrium, we reproduce the spectrum usually provided by the 
fluctuation-dissipation relation. Fluctuations in unstable systems are also discussed and, as an example, we compute 
the fluctuation spectrum of longitudinal fleld in the two-stream system. 

The method under consideration, although physically appealing, is certainly not the most effective to analyse 
equilibrium plasmas. And our actual goal is to set a stage for nonequilibrium calculations similar to those of the two- 
stream system. Our particular interests is focused on the quark-gluon plasma - a highly relativistic system governed 
by nonAbelian dynamics which, in spite of important differences, manifests profound similarities to electromagnetic 
plasmas discussed at length in [8| . The quark-gluon plasma produced in relativistic heavy- ion collisions is presumably 
unstable to chromomagnetic modes, see the review [9]. The instability growth is associated with generation of 
chromomagnetic fields which in turn strongly infiuence transport properties of the plasma [lO] . The fiuctuation 
spectrum of chromomagnetic fields is an important issue to be settled. 

Our paper is organized as follows. In Sec. |TT] we present the theoretical framework to be used in our further 
considerations. The linearized kinetic equation are solved together with Maxwell equations by means of the one-sided 
Fourier transformation in Sec. IIIII The electric and magnetic fields are expressed through the initial values of the 
fields and electron distribution function. Sec. IIVI deals with the initial fiuctuations. Those of the distribution function 
are identified with the fiuctuations in a classical system of noninteracting particles. The initial fiuctuations of fields 
are expressed through the particle fluctuations using the Maxwell equations. The well-known fluctuation spectrum 
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of longitudinal electric field is obtained in Sec. |V] while the fluctuation spectra of magnetic and electric flelds are 
derived in Sees. IVII and lVIII respectively. It then becomes clear why the analysis of longitudinal electric field is much 
easier than that of the general case. In Sec. IVIIII we extend our calculations to nonequilibrium anisotropic plasma, 
discussing fluctuations of longitudinal electric fleld in the unstable two-stream system. Our results are summarized 
and concluded in SecHX] Throughout the article we use the CGS natural units with c = fcs = 1. 



II. PRELIMINARIES 



We consider a classical plasma where ions are assumed to be a passive background which merely compensate the 
charge of electrons. However, the ions can be easily included in the considerations. The time scale of fluctuations of 
interest is much shorter than that of inter-particle collisions, and consequently the starting point of our analysis is 
the collisionless transport equation of electrons 

((^ + V • V) - e(E(i, r) -f V X B(t, r)) • V,,)) /{t, r, p) = , (1) 

where /(t, r, p) is the distribution function; E(i,r) and B(t,r) denote the electric and magnetic fields in the plasma. 
The transport equation ([T|) is supplemented by the Maxwell equations 

V-E(t,r) = 47rp(t,r), V-B(t,r)=0, (2) 



dt ' ' ' dt 

with the electric charge density and current given as 

(Pp 



p(t,r) = -e I j^^f{t,r,p) +enionB , (3) 



= -e y"^|-^v/(t,r,p) 



(2^)3-' ■ J (2^)3 

The transport equation linearized in 5/ is 



and 



f d'^p 

p(t,r) = j^5f{t,r,p) , 

d^p 

W) 



(4) 



The distribution function is assumed to be of the form 

/(t,r,p) = /0(p)+<5/(t,r,p) , (5) 

with 

/°(p)»<5/(t,r,p), |Vp/"(p)| » |Vp<5/(i,r,p)| , (6) 

and 



(7) 



( + V • V)5f{t, r, p) - e(E(t, r) + v x B(i, r)) • Vp/" (p) = . (8) 



(9) 



^ -ej -0^^Sfit,r,p) . (10) 
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III. INITIAL VALUE PROBLEM 

We are going to solve the linearized transport equation ([8]) and Maxwell equations ^ with the initial conditions 
5/(t = 0,r,p) = d7o(r,p) , E(t = 0,r) =Eo(r) , B(i = 0, r) = Bo(r) . (11) 
We apply to the equations the one-sided Fourier transformation defined as 

fiio, k) = l°°dtj dVe*(-*-'' "-)/(<, r) . (12) 

The inverse transformation is 

where the real parameter cr > is chosen is such a way that the integral over u is taken along a straight line in the 
complex w— plane, parallel to the real axis, above all singularities of /(a;,k). 
We note that 

' dt I dVe'(-*-''--) = ^iu;f{u;, k) - f{t = 0, k) . (14) 

The linearized transport ([8]) and Maxwell equations ([2]), which are transformed by means of the one-sided Fourier 
transformation, read 

-i{Lo~k- ^^)^fiu, k, p) - e(E(c^, k) + V X Biu, k)) • Vp/°(p) = 5/o(k, p) , (15) 

ik • E(t^,k) = 47rp(w,k) , ik-B(w,k) = 0, (16) 

ik X E{oj, k) = iujB{uj, k) + Bo(k) , ik x B(w, k) = 47rj(w, k) - icjE(tj, k) - Eo(k) . 
One finds the solution of the transport equation as 

<5/(u;, k, p) = {e{E{Lu, k) -f v x B(c., k)) • Vp/"(p) + 6fo{k, p)) . (17) 

A. Electric field 

Substituting the solution ^T7\ into the Fourier transformed current pH)) and using the third Maxwell equation 
to express the magnetic field through the electric one, the current gets the form 

j(^,k) = -^e^J ^__^((l-il-I)E(^,k) + l(v.E(^,k))k).V,Ap) (18) 

+ ^'/7?T3^^(--xBo(k)).V,/°(p)-^e/^^^<5/o(k,p). (19) 

Since the dielectric tensor e'-' (cj,k) in the coUisionlcss limit equals [l| 

47re2 f (fp k-w^:f, v^k^\^ ^ 



the current can be written as 



r(Lo,k) = ^(e'J(a;,k)-<5^-'")i?^(c^,k) (21) 
+ /7fTJ^^-((-vxBo(k))V^/"(p)-^e /^^^5/o(k,p). 
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Combing the third and fourth Maxwell equations (|16p . one finds 

[{uj^ ~ k"^) S'^ + k'k^E^ {Lu,k) = -ATTiiu f{uj,k) + iujE'i^Ck) - i{k x Bo(k))' 
Substituting the current (|2ip into Eq. (P^ . one obtains 



-Aneoj 



(27r)3 w - k • V 
Denoting the matrix in left-hand-side of Eq. (PH)) as 

E^^ (w,k) = -k^J'J' -I- fc^F -f w2£*^'(a;,k) 
the electric field given by Eq. can be written down as 

d^p (S-i)'J(w,k)wJ' 



(5/o(k,p)+ic^£;^(k)-z(kxBo(k))' . 



(22) 



(23) 



E\uj,k) = - 47re 



vxBo(k)) • Vp/°(p)+c^5/o(k,p) 



(24) 



(25) 



(27r)3 oj-vk 

+ iiu{i:-^y^{Lj,k)Ei{k)^i{j:-'y^{u;,k){kxBo{k)y , 

which is the main result of this section. 

When the plasma stationary state described by /"(p) is isotropic, the dielectric tensor can be expressed through 
its longitudinal and transverse components 



£^^"(c.,k) =ei(c^,k) — +eT(c.,k) (S^ 
where £^(0;, k) and £t('-^, k) are well known to be equal to 



£L(w,k) = 1 



47re^ [ d^p 



dfip) 



k2 J (27r)3 w - k • V + iO+ dp ' 



£t(w, k) = 1 + 



27re2 f d^p 



1 



Lu J (27r)3 - k • V + iO+ 
The matrix S*^ (ijJ, k), which then equals 



dfjp) _ k^Vj^ ^ df{p) 



dp 



k2 



9p 



S*^'(cc;,k) =a;2£L(a;,k)- 



fc*fcJ 

k2~ 



^£T(c^,k)-k2)^5^^"- 



can be inverted as 



(S-i)^^(u;,k) = 



1 



1 



CL'2e^(ti;j k) k^ uP-£T(yi^k}j — k? 



(26) 



(27) 



(28) 



(29) 



(30) 



When the momentum distribution /°(p) is isotropic, Vj,/°(p) ^ p, and consequently (v x Bo(k)) • Vp/°(p) = 0. 
Therefore, the first term in the right-hand-side of Eq. (P5|) vanishes and the electric field is found as 



£'*(w, k) = - 47reLJ 



1 fc'fcJ 



1 



lcj2e^(a;^k) k^ iJ^ETiy^^—'^ 
1 fc^fcJ 1 



(27r)3 tj-k- V 



<5/o(k,p) (31) 



I cj2£^(ti), k) k^ a;2£7-(ti;, k) — k^ 
If the field is purely longitudinal, 



E(u;, k) = (k • E(^, k)) ^ , Eo(k) = (k • Eo(k)) ^ 



k2 ' 
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Eq. (I3T|) gives 



, ,N 47re f (Pp k-v , ik • Enfk) , , 

k.E^,k = U-f- ^5/ok,p+ 32 



Taking into account that 

zk • Eo(k) = 47rpo(k) = -4^e j ^ <5/o(k, p) , 

Eq. (15^ can be rewritten as 



1 w( ^^ 47re / d^p ^/o(k,p 

k-Ew,k= ^ — U -PT^ i • 33) 

Eq. can be obtained directly by substituting the solution of transport equation P7)) (with B = 0) into the first 
Maxwell equation. Then, the initial electric field does not show up. 

B. Magnetic field 

Using again the third Maxwell equation (fT6|) to express the magnetic field through the electric one, Eq. ((25|) 
immediately provides 



(34) 



B\^,k) ^ -e^^''k^{j:-^r{u^,-k)( - 4me' [ ^' (v x Bo(k)) • V,/°(p) 

-47rec. y ^-^L_SfoiKp)+zu;El,{k)~z{k x Bo(k))'^ + -^S^(k) . 
When the plasma stationary state is isotropic and (E^^)*^ (tj, k) is given by Eq. (|5D|). one finds 

fijl hi 

-"''^^(^")"(-''-)- .%.(.,k)-k^ - 

The first term in the right-hand-side of Eq. (f34|) vanishes, because (v x Bo(k)) • Vj,/'^(p) = 0, and thus 

^/ ,N 47re /■ (fp kxv ^ „ , ik x Eo(k) iw£T(aj,k) ^ , 

BO^.k = - , , ■ , ■ , / — ^ -— ^/o k,p + , , °y,2 + 2 / V v2 BoW- 36 

a;^£T(w, k) — k^ 7 (ztt)'^ cj — k • V cj-^et^cj, k) — k^ cj^eri'jJ, k) — k^ 

IV. INITIAL FLUCTUATIONS 

The correlation functions of electric or magnetic fields, (£''(ti,ri)i?^(t2,i"2)}, {B^{ti,ri)B^ {t2,r2)) ((■■■) denotes 
averaging over statistical ensemble), are determined by the fields E{t, r), T3{t, r) found in the previous section and the 
initial correlations (5/o(ri, pi)5/o(r2, pa)), (El{ri)E'„ir2)), (BS(ri)B^(r2)), ((5/o(ri, pi)£;i^(r2)), ((5/o(ri, pi)B^(r2)), 
and {E'^{yi)Bq{y2)) which are discussed in this section. 

We identify the initial correlation function ((5/o(ri, pi)(5/o(ri, pi)) with the correlation function 
((5/(ti, ri, pi)(5/(t2, r2, P2))free taken at ti = f2 = of the system of free classical particles (obeying Boltz- 
mann statistics) in a stationary homogeneous state described by the distribution function /*'(p). As well known 

0, 

(^/(il, n, pi)(5/(t2, r2, P2))frcc = (2^)35(3) _ p^) J(3) (j.^ _ n - Vi(t2 - h)) /"(pi) • (37) 

Then, 

(5/o(ri,pi)(5/o(ri,pi)) = {5f{h = 0, n, pi)(5/(t2 = 0, r2, P2))frcc = (2^)^5(3) (p^ _ p,) S(^\v2 - r,) /°(pi) , (38) 
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and 



(5/o(ki,pi)5/o(k2,P2)) = (27r)35(3)(p2_p^)(27r)35(3)(k2+ki)/0(pi) 



(39) 



The correlation function /c2((5/o(ki, pi)i?Q(k2)) can be also expressed through ((5/o(ki, pi)(5/o(k2, P2))- Using the 
first Maxwell equation, one finds 



And finally, 

fcifc^(i?S(ki)£;^(k2)) = - 



-47ri((5/o(ki,pi)po(k2)) 



(40) 



/Jo 
^ (<5/o(ki,pi)(5/o(k2,P2)) = 47rze(27r)3^(3)(k2 +ki) /"(pi). 



-16^^(po(ki)po(k2)) 



(2^) 



(27r)3 (27r)3 
3 /°(P) • 



((5/o(ki,pi)(5/o(k2,p2)) 



(41) 



When the electric field is not purely longitudinal, the computation of the initial correlations (£'Q(ri)i?Q (r2)), 
((5/o(ri,pi)i?Q(r2)) is more complicated, as the electric field Eo(r) is not fully determined by (5/o(r,p) but (5/(i,r,p) 
enters here. To compute {EI{vi)EI{y2)) , (5/o(ri, pi)£'^(r2)) as weU as (BJ(ri)B^(r2)), ((5/o(ri, pi)B^(r2)), and 
(i?o(ri)-Bo(r2)), we use the Maxwell equations transformed using the Fourier transformation not the one-sided Fourier 
transformation. Actually, the Fourier transformed Maxwell equations are very similar to the one-sided Fourier trans- 
formed Maxwell equations (|16p . The initial electric and magnetic fields are simply absent in the former ones. However, 
it should be clearly stated that the one-sided Fourier transformation is not mixed up with the Fourier transformation. 
The latter is used to compute only the initial fluctuations which are independent of lo. 

Combing the third and the fourth Maxwell equation, one gets the equation as Eq. ([^5)1 but the terms with Eo(k) 
and Bo(k) are absent. Inverting the matrix in the right- hand-side of the equation, we get the electric field expressed 
through the current 



i?^(^,k) 
The magnetic field is given as 



-47riw 



1 PfcJ 



j2 -k2 



f (^,k) 



B(^,k) = - 



Am 



k X j(cj,k) 



The correlation function (i?Q(ki)£'Q(k2)) is derived as 

duji du}2 



{El,ik,)Ei{k2)) 



27r 27r 

1 klk'l 
uji k? 



{E\uiM)E\L02,k2)) = -(4^)^ 



diVi duj2 



(42) 



(43) 



(44) 



1 P i-' 

UJ2 kl 



L02 / , _ fc|M \ 



k^ 



X (/(wi,ki)j^(tJ2,k2)) 

2 f duJi (hj2 (Ppi (Pp2 k 



2 



uji kj oj 



k? 



1 P 

1 ^21^2 



+ 2 .2 ^ -42^ f2^ + ^ -f2^ 

- k^ V kf / UJ2 ki - ko V kf / 



k? 



X ('5/('^i,ki,pi)(5/(a)2,k2,P2)) . 
As previously, we identify ((5/(u;i, ki, pi)(5/(tj2, k2, P2)) with (J/(wi, ki, pi)5/(w2, k2, P2))frco which equals 

((5/(c^i,ki,pi),5/(a;2,k2,P2))frcc = (27r)35(3)(p2_pi)27r(5(c^i+c^2)(27r)3,5(3)(ki+k2) (45) 

'LUi -UJ2 ki - k2 



X 2n6 



vi r(pi) 
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Then, after performing trivial integrations, (_BQ(ki)£'Q(k2)) equals 

r(p) 



i2nr' ((ki.v)2-k?)((k2-v)2-k2) • 

Computing fc5^fc2(-Bo(ki)i?o(k2)}, one reproduces the result (j41]l . Analogously to the correlation fimction 
(i;S(ki)£;^(k2)), one finds 

(ki • V2)?;2 - k\ 



(£;^(ki)5/o(k2,P2))=4^*e (27r)35(3)(j^^_^j^^)^0(p^) 



(ki • V2)2 - k? 



Starting with Eq. (|43|) . we obtain 

(i3j(ki)i3^(k2)) = -(47re)^(27r)35(3)(k^ +k2) e»'='e^™"fcfefc2" / /°(P) 



and 



Finally, one computes 



(27r)3 ((ki.v)2-k?)((k2-v)2-k2; 



(i3^(ki)(5/o(k2,P2)) =47rze(27r)3j(3)(ki+k2)r(p2) ■ 

(ki • V2)^ - kf 



(i.S(k.)i?^(k2))^-(4.ef(2.)3.(3)(k,+k2) / ^/°(P) (('^--K-^O^^-^'^.V 



(2^)3^ ((ki.v)2-k2)((k2-v)2-k 



(47) 



, (48) 



(49) 



(50) 



V. FLUCTUATIONS OF LONGITUDINAL ELECTRIC FIELD IN ISOTROPIC PLASMA 



We first consider a special case of purely longitudinal electric field in the isotropic plasma when the electric field is 
jiven by Eq. Then, 



k\ki{E\uj^M)E^{u^2M)) 



1 



uiiUJ2eL (t^i , ki ) El (cj2 , k2 ) 



(51) 



1« 2 2 / d^Pl '^^P2 ki-Vi k2-V2 

X 167r^e^ / 7:^^7:^^ ^ ^ (5/o(ki, pi)5/o(k2, P2)) 



Ame 



(27r)3 (27r)3 tji — ki • vi 1^2 ^ k2 • V2 
d3pi ki • vi 



(27r)3 LJi — ki • vi 

d3p2 k2 • V2 



(27r)3 u;2 - k2 ■ V2 

Substituting the formulas of initial fluctuations (|39l l40l l41]l into Eq. ([5T|l , one finds 



A:^(<5/o(ki,pi)i?^(k2)) 

A:l(i?S(ki),5/o(k2,p2)) - fclfc^(i?^(ki)i?^(k2)) 



A:lfc^(S^(wi,ki)£;^(c^2,k2)> - 1671^6 



2„2 



(2^)3,5(3) (k2+ki) 



Wl W2eL (wi , ki ) £l (W2 , k2 ) 



(52) 



(2^ 



ki ■ V k2 • V ki • V 



ko • V 



Wl — ki • V 0^2 — k2 • V CJi — ki • V W2 — k2 ■ V 



/°(p) . 



It appears that 



ki • V k2 • V 



ki • V , k2 • V _^ _ tJi 



CJ2 



0^1 — ki • V C1J2 — k2 • V — ki • V 0^2 — k2 • V lji — ki • v 0^2 — k2 • v 

and consequently Eq. ([52ll simplifies to 



fcifc^(i;*(a;i,ki)i;J(w2,k2)) = 16^ 



2„2 



(2^)3^(3) (k2+ki) /■ <Pp /0(p) 



eL(t^i, ki) eL(t^2, k2) 7 (27r)3 {ui - ki • v)(w2 - k2 • v) 



(53) 
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This equation could be easily obtained directly from Eq. ([33]) where the initial electric field is already eliminated. 

Taking into account that the electric fields are parallel to their wave vectors, and consequently 
kiki{E^{uJi,ki)E^ {102,-^1)) = kl{E'{uJi,ki)E'{uj2, -ki)), one finally finds 



2 2 (2^)3^(3) (k2+ki) f <fp r(p) 



kfez,(tJi, ki) £^(^2, -ki) J (27r)3 [uji - ki • v){w2 + ki • v) 
Let us now compute ri)£"(i2, i"2)) given by 

Zeros of ei,(wi,ki) and of the denominators {iHi — k^ • v + iO+) with j = 1,2 contribute to the integrals over uji and 
tj2- However, once the plasma system under consideration is stable with respect to longitudinal modes, all zeros of 
El lie in the lower half-plane of complex lj. Consequently, the contributions associated with these zeros exponentially 
decay in time and they vanish in the long time limit of both ti and t2. 

We are further interested in the long time limit of (£''(ti, ri)i?*(t2, r2)) and then, the only non- vanishing contribution 
corresponds to the poles at wi = ki • v and 1^2 = k2 • v. This contribution reads 

m.M.^.f,^ = '^^'1 ^ I (l^--"''''--'-'-'"-"'> .V.(..v,ff(-.<.v.-K) ■ 

Keeping in mind that £l(— cj, — k) = s*j^{uj, k) for real to and k, it can be rewritten as 

As seen, ri)i?'(i2, r2))oo given by Eq. ([57]) depends on ti, t2 and ri, r2 only through {ti — 12) and (ri — r2) and 

it can be written as 

{E\t,,r,)E\t2.V2))oo = f '^^e-<-('^-^-)-^<^^-^-)){E^E^)^^ , (58) 
where the fluctuation spectrum {E'^E^)u]k is 

= k^M^ y - • • 

In the case of equilibrium plasma, the formula (|59p provides the result which can be obtained directly by means of 
the fluctuation-dissipation theorem. Let us derive the result. Due to the identity 

^ r- ip iTTS{x) , (60) 



the imaginary part of SLiujjk), which is given by Eq. (|27p. equals 



^eU^, k) = y ^ ^(. - k . V) k . . (61) 



In equilibrium /'^(p) e ''■^p and 3/*'(p)/9p = — /9v/°(p). Therefore, Se^ equals 

47r^e^ f d^p \ j-or \ 47r^e^ti; /" d^p 



(27r)3 ^ ' ^ ' rk2 y (27r)3 

Using the expression ([S^ . the formula (|59p is rewritten as 



(£;'i?':.k = 8vr^^^f^, (63) 
|eL(w,k)|2 



which agrees with Eq. (51.25) from [7] which is obtained there in essentially the same way. 
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VI. FLUCTUATIONS OF MAGNETIC FIELD IN ISOTROPIC PLASMA 



As seen in Eq. (|36p . the magnetic field in isotropic plasma is given by three terms. Therefore, nine terms enter the 
correlation function (i?'(a;i, ki)i3-' (^2, k2)}. Substituting into these terms the initial fluctuations derived in Sec. llVi 
one finds after an elementary but lengthy and tedious analysis the following expression 



{colsTiioM - kf )) {colsT{c02,k2) ~ k^)) 



(64) 



(27r)3 ' {cj, - ki • V)(C.2 - k2 • V)((ki • V)2 - k2)((k2 • v)^ - k^) 

X (wi(ki ■ v) - kj) + wieT(wi,ki)(wi - ki • v) 

X (cJ2(k2 • v) - k2) + UJ2eTiuJ2, k2)(w2 - k2 • v) 

We now compute ri)i3^(<2, i"2)) given by 



/m/, ^m7. ^\ r^'^ dLo^ r^'^ duj2 [ Sk^ [ <Pk2 



(65) 



X e 



Zeros of {lofeTi^i, k^) — k|) and of oj^ — k^ • v + iO^) with i = 1, 2 contribute to the integrals over uji and lu2- However, 
once the plasma system under consideration is stable with respect to transverse modes, all zeros of (ujfeTi^^i, ki) — k|) 
lie in the lower half-plane of complex uj. Consequently, the contributions associated with these zeros exponentially 
decay in time and they vanish in the long time limit of both ti and t2- 

We further consider the long time limit of (i?*(ti, ri)i?^ (^2, r2)) and then, the only non-vanishing contribution 
corresponds to the poles at = ki • v and UJ2 — \<.2 ■ v. This contribution reads 



{B{t„v,)B {t2,V2))^ - ~J ^^^f (p)e 

(47re)2(27r)3,5(3)(ki -f k2) e*'='eJ™"A:f fc^ 



(66) 



((ki.v)2-k?)((k2-v)2-ki) 

It can be easily expressed as 

{B\h,r,)B'{t2,r2)) 
where the fluctuation spectrum is 
{BW^)^i, 



IsTicJ.M) - k?)) {cj^,eTiu;2,k2) - k^)) 

(coi(ki.v)-k2)(c^2(k2-v)-k2) 



duj (fk 
2^(27r)3' 



327r3e2e*"e-'™"fc'=fc™ 



i(i^(ti-t2)-k-(ri-r2)) (^B^B^)^i^ 



"^^^ fip) S{iJ-k- v) w'w" 



{cu'^eTiuj, k) - k2)) {io^eri-uj, -k) - k2)) J {2tt)^ 
When both lu and k are real ETi—i-^, — k) = s^{u!, k). Therefore, the fluctuation spectrum can be rewritten as 

327r3e2e^'=^'eJ™"/c'=fc" f (fp ^ 



(B'B^)^, 

liij^eTi^^, k) — k2 

One observes that the matrix function 



(2^)= 



■ f{p)d{uj -k-v) v'v 



(2^) 



(67) 



(68) 



(69) 



(70) 
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which enters the correlation function ((69)) . can be decomposed as 



M^' [u, k) = MLiuj, k) -j^ + A/t(co, k) (^'^ - — j , 
because the plasma is assumed to be isotropic. Comparing Eq. (|7T|) to Eq. (|7D|) . one finds 



Mi(w,k) = 



3 /"(P) <5(w - k • v) 



(2^) 



k2 



Mt(c., k) ^ - / ^ /0(p) - k . v) 



(k.v)2 

2 7 (27r)3 

Using the decomposition (|7ip . the correlation function (|69p can be written down as 



Iw^erCw, k) — k^l 



-Mt(w, k) 



(71) 



(72) 



(73) 



(74) 



For equilibrium plasma the correlation function {B^B^)i^\^ can be expressed in the form of fluctuation-dissipation 
relation. One first observes that due to the identity ([SO]) , the imaginary part of £t('j-', k), which is given by Eq. 
is 



9£t(w, k) 



LO J (27r)3 

With the equilibrium distribution function, SJet equals 

27r2e2 f (Pp 



5{u) — k • v) 



■ 9/°(p) k-v ar(p) 

V • — -IT- k 



9p 



k2 



9p 



Twk2 7 (27r)3 

Consequently, the function Mt ([73]) can be expressed through Set ([76|) as 



Tuj 
47r2e2 



3eT(w,k) , 



and finally, 



eT(w,k) 



k2 |2 



(75) 



(76) 



(77) 



(78) 



Eq. ([78]) coincides with the formula (11.2.2.7) from obtained there directly from the fluctuation-dissipation theorem. 



VII. FLUCTUATIONS OF ELECTRIC FIELD IN ISOTROPIC PLASMA 

The analysis of electric field fluctuations is much more complicated than that of the magnetic field. First of all, there 
are five terms which enter the formula of electric field given by Eq. (pij) . and consequently, the correlation function 
(i?*(ci;i, ki)i?^ (aj2, k2)) includes 25 terms. The magnetic field is purely transverse and some terms automatically drop 
out but the electric fields have longitudinal and transverse components. Using the formulas of initial fiuctuations, 
which are derived in Sec. IIV| and patiently analyzing term by term, one obtains after an elementary but very lengthy 
calculation the correlation function of the form 



{E\lo^M)E'{^2M)) 



167r2e2(2.)3p)(ki+k2) j ^ f{p) 



(79) 



k\ 



jfeL{i^i,'ki) w|eL(a)2,k2) kf (wi - ki • v) k^(w2 - k2 • v) 



kl 



v^kl - k{{\<.2 ■ v) UjI[uJ2{uJ2{^2 ' v) - k|) - k|(w2 - k2 • v) 



£l(wi, kl) w|eT(t^2, k2) - k| kj(wi - ki • v) k^(w2 - k2 ■ v) ((k2 • v)2 - kg) 
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1.3 

^2 



a;|[ijji(a;i(ki • v) — k^) — k^([Ji — ki • v) 



ujferiuji^'iii) - k^ wfg^^^^, k2) k2(wi - ki • v) ((ki ■ v)2 - kf) kl{uj2 - k2 • v) 
k\ (ki • v) - v^kj kl (k2 • v) - vi^l 



W^£T(wi,ki) - kf w|£T('^2,k2) - k| 

cji(a;i(ki • v) k^) - kf (tJi - ki • v) cj2(w2(k2 • v) - kj) - k|(tj2 - k2 • v) 



k2(c^l - ki • V) ((ki • V)2 - k?) 



We now compute ri)£'-' (^2, i"2)) given by 



(i?^(ti,ri)£;^(t2,r2)) = ^ ^ 



27r 



27r 



k2(c^2-k2-v)((k2-v)2-k2) 



(Pki f (Pk2 



(27r)3 7 (27r)3 



(80) 



^-i(ii;iti-ki-ri+ct;2t2-k2-r2) 



(£;'(^i,ki)i;j(w2,k2)) 



Zeros of (wf eT(wi, ki) — k^), cj|£i(a;i, k^) and of [uji — k^ • v + iO+) with i — 1,2 contribute to the integrals over uji and 
UJ2- However, once the plasma system under consideration is stable, all zeros of (a;feT(wi,ki) — kf) and w?eL(wi,ki) 
lie in the lower half-plane of complex u. Consequently, the contributions associated with these zeros exponentially 
decay in time and they vanish in the long time limit of both ti and ^2 ■ 

We further consider the long time limit of {E^{ti,ri)E^ (12,^^2)) and then, the only non-vanishing contribution 
corresponds to the poles at wi = ki • v and u;2 = k2 • v. This contribution reads 



{E^h,n)E^it2,r2))c 



(fki (fk2 

(2^(2^ 



(81) 



fO( \ -i(^iti~ki-ri+i^2t2~k2-r2) ^l'-^2 

(27r)3 ■' km 



ujik\ 



fc*(ki • v) - v'k? 



w^£L(wi,ki) wJeT(wi,ki) - k^ 
The correlation function (ISTI) can be rewritten as 



U!2kl 



^ ^^•kl-fc^(k2-v) 

wi£L(w2,k2) w|£T(w2,k2) - k| 



cJi— ki-v, , a;2— k2-v 



{E\tuYi)E^{t2,Y2))o. 

where the fluctuation spectrum is 

{E'E^)^^, = l(!>^''e^ 
{ k' k3 



duj d?k 
2^(2^' 



-i(w(ti-t2)-k-(ri-r2)) l^J^ij^j"^ 



UJ k ; 



(27r)3 



fip) 27rJ(c.-k.v) — 



k' 



k4 

v^k"^ -y(k-v) 



(82) 



(83) 



+ 



w2e^((jj,k) w2£^(— oj^k^) 
v'k"^ - k'{k-v) k'-i 



ct;2eL(aj,k) a;2eT(— w, k) — k^ 

fc^(k-v)-w'k2 F(k-V) - W^k^ 



j2eT(w,k)— k^ cj2£^(— k) 



;2eT(^,k)— k^ cj^g^^^^^^ j^^ _ jj2 



One easily proves that the second and third contribution to the fluctuation spectrum (|83p vanish due to the plasma 
isotropy. Taking into account that for real uj and k, the dielectric function obeys es(— cj, — k) = s*{lo, k) with s = L, T, 
the fluctuation spectrum (j83p can be written as 



(E'E^)^), - 16^'e' 



(P(k. 



^^^^3Ap)2..(.-k.v)^ 



(84) 



'*k2)(fcJ(k- V) - W^k^) 



j2£T(w,k) -k2|2 
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Due to the plasma isotropy, the expression, which enters the transverse contribution, can be further rewritten as 



/°(p) 2tt5{uj - k • v)(fc'(k • v) - v'\^) (F (k • v) - w^k^) 



(2^) 



k2 \uj^eL{^,k)\ 



3eT(t^,k) 



(85) 



In the equihbrium plasma, the imaginary parts of eL(a-', k) and exi'-^^k) are given by the formulas (|^ [75]) and the 
fluctuation spectrum (|84p can be expressed as 



(86) 



which for the longitudinal fields reproduces the formula ([55)) . The result agrees with Eq. (11.2.2.6) from [H 
derived using the fluctuation-dissipation relation. 



VIII. FLUCTUATIONS OF LONGITUDINAL ELECTRIC FIELD IN THE TWO-STREAM SYSTEM 

Nonequlibrium calculations are usually much more difficult than the equilibrium ones. The first problem is to invert 
the matrix Yi^^{Lo, k) defined by Eq. ([24ll . In the case of longitudinal electric field, which is discussed here, it is solved 
trivially. We start with Eq. (|23p projecting it on k and assuming that E and Eq are purely longitudinal. Then, the 
matrix ([24)1 is replaced by the scalar function. 

Further, we neglect the first term in the r.h.s. of Eq. ([23]) . This term vanishes in isotropic systems; it is of order 
higher than the second term; it is also expected to be small in nonrelativistic regime due to the smallness of particle 
velocity. So, there are good reasons to neglect it. Eliminating Eq by means of the first Maxwell equation we obtain 
Eq. ([33]) which was previously derived for the case of isotropic plasma. In the following we consider fluctuations of 
longitudinal electric fields in the nonrelativistic two-stream system. 

The distribution function of the two-stream system is chosen to be 

/"(p) = (2^)3n[<5(3)(p_q)+^(3)(p + q)] ^ (gy) 

where n is the electron density in a single stream. To compute ei(a;,k) we first perform integration by parts in 
Eq. (j27|) and then, substituting the distribution function (|87|) into the resulting formula, we obtain in the nonrelativistic 
approximation 

{u-^{k■nfY-2^?{^- + {k■uf) 

£L(,t^,K) — —2 , 



(a;2 - (k ■ u) 



2\ 



where u is the stream velocity (nonrelativistically u — q/m with m being the electron mass) and /i^ = Ane^n/m. 
There are four roots ±ct'±(k) of the dispersion equation ei(w, k) = which read 



Lol{k)^ + {k-uf ±^i^j:flTW^ ■ (89) 

As seen, < Ci-'+(k) £ R for any k but ijJ-(k) is imaginary for (k • u)^ < 2/^^ when it represents the well-known 
two-stream electrostatic instability. For (k • u)^ > 2/1^, the mode is stable, < t^J-ik) G -R- 
The correlation function (i?*(a;i, ki)i?*(a;2, k2)) as given by Eq. (|54p equals 

{E\ujiM)E\uj2, k2)) = -IQ^^e'n^^^^^^^-l^^^^ [(c^i -)- ki • u){uj2 -)- k2 • u) -f {oji - ki • u)(u;2 - ka • u) 

- (ki • u)2 



(cji - c^-(ki)) (wi -1- ti^-(ki)) (wi - w+(ki)) (tJi + ti^+(ki)) 

- (k2 ■ u)2 ^ 

(cj2 -cj-(k2))(a;2 + t^-(k2)) (tJ2 -w+(k2))(a;2 -l-a;+(k2)) 



(90) 
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One observes that the poles of the correlation function (i?*(cL'i, ki)i?'(w2, k2)) at wi = kiv and uj2 = k2V, which 
give the stationary contribution to the equilibrium fluctuation spectrum, have disappeared in Eq. (j90p as the inverse 
dielectric functions vanish at these points. 

The correlation function {E^{ti,ri)E'{t2,r2)) is given by Eq. ^ with ki)£;'(cj2, ka)) defined by Eq. H^ . 

Performing the trivial integration over k.2 and taking into account that a;±(— k) — a;-|-(k), one finds 



pOO-\-i(7 


did I 


noo-\-ia 


duj2 f 


' — oo-\-ia 


27Ti ^ 




J 



r)3 



k2 



■ [wiW2 - (k • u)2] 



(k-u)2 



[uJi - t^-(k)) [uji + w_(k)) (wi - w+(k)) (uji + uj+ik)) 



(k-u)2 



[uj2 — t^-(k)) {uJ2 + ^-(k)) (cLi2 — t^+(k)) (0^2 + t^+Ck)) 



(91) 



There are 16 contributions to the integrals over lox and L02 in Eq. (j9ip related to the poles at ±aj±. Summing up the 
contributions, we get after lengthy calculation 



{E\ti,vi)E\t2,V2)) = 16^Vn 
K-(k-u)2)^ 



1 



(k.u)^)(c^^-(k-u)^ 



(cj2 - (k • U) 



(27r)3 k2 (w2_^2)2 

(tj![ - (k • u)2) cos {uj+{ti + t2)) + {cjI + (k • u)2) cos {uj+ih - 12)) 

(bj^bj^ — (k ■ u)^) cos(ti)+ii + uj^t2) + + (k • u)^) cos(a;+ti — a;_t2) 

+ {lj^u}^ — (k ■ u)^) cos(cj_ti + uj+t2) + + (k • u)^) cos(a;_ii — cj+t2) 

{ul - (k • uf) cos + t2)) + (tJ^ + (k ■ uf) cos - t2)) 



(92) 



Let us now consider the domain of wave vectors (k ■ u)^ < 2/^^ when w-(k) is imaginary and it represents the 
unstable electrostatic mode. We write down a;_(k) as with < 7k £ i?, 



7k 



tl^/Jj?+I^~Vi^ - - (k • u)2 



(93) 



We are interested in the contributions to the correlation function coming from the unstable modes. The contributions, 
which are the fastest growing functions of (ti + ^2) and (^i — 12), correspond to the last term in Eq. (j^ . The 
contributions provide 



unstable 



47r^e^n 



pik(ri-r2) 



(7^ + (k-u) 



(27r)3 k2/i2(Ai^+4(k-u)2) 
(7^ + (k • uf) cosh (7k(ti + t2)) + (7k - (k • u)2) cosh (7k(ti - ^2)) 



(94) 



where we have taken into account that uj\_ — lo seen, the correlation function (|94p is 

invariant with respect to space translations ~ it depends on the diff'erence (ri — r2) only. The initial plasma state is 
on average homogeneous and it remains like this in course of the system's temporal evolution. The time dependence 
of the correlation function (I94p is very different from the space dependence. The electric fields exponentially grow 
and so does the correlation function both in (ti + 12) and (ti — t2). The fluctuation spectrum also evolves in time as 
the growth rate of unstable modes is wave- vector dependent. After sufficiently long times the fluctuation spectrum 
is dominated by the fastest growing modes. It should be remembered, however, that our results hold for times which 
are not too long. Otherwise, the perturbation, which exponentially grows, violates the condition ^ justifying the 
linearization procedure. 
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IX. SUMMARY AND CONCLUSION 



The calculations presented here show how to obtain a spectrum of electromagnetic fluctuations in equilibrium or 

noncquilibrium plasmas as a solution of an initial value; problem. Wc first linearize the transport equation around 
the state which is on average stationary and homogenous. The linearized transport equation is solved together with 
the Maxwell equations by means of the one-sided Fourier transformation. The time dependent fluctuation spectrum 
is expressed through the fluctuations in the initial state. Electromagnetic initial fluctuations are determined by the 
initial fluctuations of the distribution function. The later are identified with the fluctuations in a classical system 
of noninteracting particles. We compute fluctuation spectrum of longitudinal electric flelds in isotropic plasma, and 
then there are considered fluctuations of magnetic and electric fields. Our equilibrium results coincide with those 
obtained by means of the fluctuation-dissipation theorem. However, the method adopted here clearly shows how 
the system looses its memory and how the stationary equilibrium spectrum of fluctuations emerges. As an example 
of unstable systems, the fluctuations of longitudinal electric fleld in the two-stream system are considered. The 
fluctuation spectrum appears to be qualitatively diff'erent than that of the equilibrium plasma. 

The scheme of calculation, which is worked our here in detail, can be applied to a variety of plasma noncquilibrium 
states. Our actual objective is, however, to generalize the scheme to study fluctuations in the quark-gluon plasma 
mentioned in the Introduction. Such a generalization is not quite trivial. The problem is that chromodynamic 
flelds are, in contrast to their electromagnetic counterparts, gauge dependent, but physically meaningful correlation 
functions have to be gauge independent. A treatment of color charges needs to be different as well. 
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